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ABSTRACT 


In this thesis an attempt is made to determine the viscous 
effects and contribution of pressure gradient terms to velocity- 
profiles in two-dimensional, turbulent, submerged jet using the 
perturbation technique* The problan. of two-dimensional jet is 
shown to be one of regular perturbation* She reciprocal of the 
Reynolds number in the vorticity transport equation for the two- 
dimensional turbulent flow is used as the perturbation parameter. 
Jrandtl * s eddy viscosity model which takes the eddy viscosity to 
be constant across the jet-width is used as the closure hypothesis. 

The stream function, Reynolds stress and viscous stress terms are 
expanded in asymptotic series. We obtain the first order and 
second order equations by comparing the orders of the various terms 
in the vorticity transport equation. 3he first order equation is 
the boundary-layer equation for zero pressure gradient and was 
solved by Goeftler (8) for the case of two-dimensional jet. The 
second order equation -vdaich includes viscous effects is solved 
numerically. The pressure gradient term is found to be of a still 
higher order and thus does not affect -the velocity distribution 
to second order. 

Due to viscosity the velocity near the jet axis reduces and 
that away from the axis increases resulting in a flattening of the 
velocity profile and a increase in the width of the jet boundary layer* 



There is also a slight dependence on "the axial distance, the viscous 

corrections reducing as we move away frcm the slit. These results 

2 4 

are valid and significant in the Reynolds number range 10 < Re < 10 

and for x-stations 8 < x/d < 25. For x/d > 25 and Re > 10^ the first 

order velocity profile (based on infinitely large Reynolds number 

and full development assumptions) can be taken as valid. For 

2 

x/d < 8 and Re <10 , the jet model breaks down and the results 
obtained are no longer valid.. 
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Constant scaling factor for -velocity. 

Constant scaling factor for length.. 

Jet frolf-widttu 

Ooostaot of proportionality for eddy riacosity ooaffioiant. 
Width of jet slit. 

Sequence of gauge functions for stream function. 

Gauge faction for eddy viscosity coefficient. 

Momentum flux in the axial direction. 

Eddy viscosity coefficient. 

Mixing length. 

Reynolds number. 

Mean axial velocity component, 
fluctuating axial velocity component. 

Mean centre line velocity. 

NoiMimensiooal mean axial yelooity component. 

mean oentre line yelooity. 

Ifon-dimensional jet exhaust velocity. 

Mean transverse velocity component, 
fluctuating transverse velocity component. 

Mon-dimensional mean transverse velocity component. 

Axial coordinate. 

Norvdimensional axial coordinate. 

Scaled monr-dimensio nal axial coord mat 



Transverse coordinate. 


Non-dimensional transverse coordinate. 

Transverse distance from axis at which the axial velocity 
is half the centre line velocity.. 

Scaled non-dimensional transverse coordinate. 

Boundary layer thickness. 

A small parameter, function of Reynolds num b er. 

Stream function. 

Density. 

A non-dimensional similarity parameter. 

A small parameter. 


Reynolds stress components. 



CHAPTER 1 

IHTROIUCTION 




INTROnJCTION 


theoretical prediction of mechanism of transfer of mass, 
momentum etc. and velocity profiles in turbulent jet flows, as 
accurately as possible, is important for optimization of jet 
engine performance while minimizing noise generation. It is also 
necessary for a better understanding of mechanian of turbulent 
diffusion in general which would be useful in designing combustion 
chambers, jet pumps, mixing tanks and disposal of waste gases. 

The majority of analytical and experimental research 
oonducted on turbulent jets has been concentrated on circular jets. 
This is mainly because of the fact that circular jets are of 
practical interest in engine exhaust configurations. Never the less 
the two dimensional geanetry also offers many opportunities for 
reaching a better understanding of the problem of free-turbulent 
flows, that is turbulent flows in absence of solid boundaries. 
(Slough plane turbulent wakes have been extensively studied by 
Townsend(23), Grant ( 9 ) etc. this is not so with plane jets. 

The need to improve our understanding of the two-dimensional 
turbulent jet has also arisen with the development of short take- 
off ard larding air crafts for which aerodynamic designers are 
considering the use of various bl own— flap concepts to provide the 
lift augmentation required for these types of aircrafts. A slot 
jet above the flaps may be used in such aircrafts. 
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In the analytical formulation of free turbulent shear flows 
the viscous effects have so far been totally neglected (20). This 
gave reasonable agreement with experimental results for high Reynolds- 
number flows (8). This may not be so when lower Reynolds-numb e r 
flow® are considered. It is also worthwhile to estimate the order 
of error introduced because of neglecting pres sure -gradient terms 
in free jets. The development of perturbation techniques as 
applicable to fluid mechanics permit an attempt to answer the above 
questions. Similar: studies have in the past been conducted by 
researchers like K.S. Yajnik (24) and Noor Afzal (12) on wall 
turbulent flows (boundary layer flows, channel flow, pipe flow etc.) 
But free shear flows had somehow been left out. In this work we 
have applied the perturbation procedure to one case of free-shear 
flows, namely, two-dimensional turbulent jet. 

One of the pioneers in the field of jet flows is H.Schlichting 
who starting from the fundamental viscous, incompressible flow 
equations obtained by numerical technique the velocity distribution 
across lamirar jets, both plane and round. Analytical solutions 
were obtained later for this case by Bickley ( 3 ) and Goldstein ( 7 ). 
The velocity distribution across turbulent jets has been investigated 
by various workers. In analytic formulation of turbulent shear flows, 
the equations of motion contain too many unknows to be solved 
completely. To obtain the results one introduces some sort of a 
closure hypothesis about the nature of turbulent stresses. Most of 
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the successful hypothesis define and use the concept of a mixing length, 
The concept of mixing length was first introduced by Prandtl in 1925 . 

He defined it as that distance in a direction transverse to the flow 
which must be covered by a lump of fluid particles travelling with 
its original mean velocity in order to make the difference between 
this velocity and the velocity at the new location equal to the 
transverse fluctuations in velocity. If the momentum is preserved 
over this length, then the arrival of this particle at the new 
location will appear as a fluctuation in the velocity at that point. 

The concept of mixing length is same what analogous with that of mean 
free path in kinetic theory of gases. Other researchers have used 
variations of the same basic idea. 

Some of the more successful closure hypothesis vised for stresses 
in free turbulent shear flows are given below : 

i) Prandtl 's old theory of momentum transfer (13) 

ii) Taylor's vorticity transfer theory (l 0 ) 

iii) Prandtls new theory of momentum transfer (15) 

iv) Reichardt 1 s molecular ana logy theory. (17) 

Prandtl's old theory assumes that the momentum mixing length 
is constant in the transverse direction, but increases linearly in 
the longitudinal direction. The first assumption uses the fact that 
there are no walls in the flow field. The mixing length decreases 
near the walls (at walls it must be zero). ToUmein (22) applied 
Prandtl’s old theory to solve the following three problems 
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a) The boundary layer of an infinite plane parallel j et 

b) Parallel jet issuing from a narrow orifice 

o) Axially symmetric jet issuing from, a small orifice# 

He was able to obtain quite satisfactory velocity fields. One 
necessary conclusion of this hypothesis is the similiarity of the 
temperature and the velocity fields; this is grossly erroneous. 

Taylor sought to correct the above by assuming that the Reynolds 
stress is determined not by transfer of momentum but by transfer of 
vorticity. He defined a vorticity mixing length as the distance over 
which vorticity of fluid particles remains constant during transfer 
from one location to another. Howarth (l 1 ) applied Taylors assumption 
to determine the velocity distribution across j ets. The results were 
identical with those of Tollmeins. However the taaperature profiles 
obtained by Taylor's theory were confirmed by experiments. 

With both these theories the velocities near the centre are 
substantially in error. To over come this difficulty Prardtl later 

argued that it is not the mixing length which remains constant across 

0 ■ 

the jet but the coefficient of turbulent momentum exchange remains 
constant. Further he proposed that this coefficient of turbulent 
viscosity or momentum exchange is proportional to the density, the 

maximum velocity difference and the jet thickness. 

Goertler ( 8 ) solved the problem of plane submerged j et using 
Prandtl's new theory and obtained for the velocity distribution 
across the jet 
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2 

= ( 1 - tanh n) where 

= local velocity 
= Velocity on jet centre line 
= Similarity parameter oy/X 
= lateral distance of the point from jet 
= Distance of cross section from the jet 
= A constant determined experimentally. 


axis 

origin 


He neglected the viscous stress completely and assumed that the 
pressure is constant everywhere. 


Reichardt assumed that diffusion process for masses of fluid 
in turbulent motion is analogous to that in molecular motion. He 
postulated that the lateral transport of momentum is proportional to 
the transverse gradient of the horizontal component of momentum and 
obtained similarity in the mean value of dynamic pressure rather than 
in velocity profiles. Velocity distribution obtained by him is given 
by . § 



-const t) 


where 


Y/X. 


Some of the earlier experimental investigations on plane jets 
were carried out by Bicknell (4 )> Dorthmann ( 6 ) and Reichardt (l 6 ) 
and later by l.J. S. Bradsbury (2) . 

Rorthnann used a rectargular slit of 3 cm x 63 cm and a discharge 
velocity U of 3500 cm/ sec. The velocity at ary point was determined 



6 


from a total head tube (not corrected for error due to turbulence). 
Porthmann only measured the velocity in the plane of symmetry 
perpendicular to the slit. 

Reichardt used a rectargul ar slit of 0.7 cm x 15 cm and a 
discharge velocity of 5000 cm/sec. In these experiments the Reynolds 
numbers specified as Re = U d/v amounted to 70000 and 23000respectively. 
Prom Porthmann’ s experiment the mixing length from the momentum 
transfer theory is 

£ » c X » 0.0165 X 

Bradsbury used a slot of 3/8" x 18" and the jet exhausted into 
a slow moving stream of air. He compared results with those of 
Townsends wake investigations. It was found that self similiarity 
in velocity profiles is established about 30 jet widths downstream 
of the nozzle. 

All the investigations reported in the literature concern them- 
selves with studying the velocity profiles for high Reynolds stress 
(viscous term negligible) and for large distances away from the orifice 
(so that complete self similarity is assumed). 

This present investigation attempts to include the effect of 
viscous stresses as perturbations to the previously obtained solutions 
valid for very high Reynolds numbers. The results obtained are 
applicable to moderate Reynolds cumbers. A weak dependence on x 
is also obtained by the perturbation procedure. 
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imoHJCTION TO PERTORBATIOff TECHNIQUES 


Introduction 

Even after a mathematical model of a physical problem has 
been constructed, obtaining exact solutions is a difficult task in 
most engineering situations. The equations encountered may be 
non-linear, coupled, of very high order or may have singularities 
which make obtai n ing analytic or even numerical solutions very 
difficult. In such situations one often resorts to some or the 
other of approximation method. One powerful approximation method 
consists of the expansion of solution in a power series in a 
parameter. This is the classical perturbation technique upon which 
much of the edifice of science rests. 

The basic premise on which perturbation technique depends 
is the hypothesis that small causes produce small effects. Thus, 
if viscous forces are snail compared to other forces in the problem, 
the solution obtained on neglecting the viscosity should differ only 
a little from the exact solution. The result can then be modified 
by including the viscous forces estimated by the derivative of the 
approximate velocities so obtained. The correction thus obtained 
is termed the second order term in the expansion of results in 
terms of the viscosity parameter. An iterative or a perturbation 
scheme is thus set up. Hie first approximation is thus valid for 
very low viscosity effects (very high Reynolds numbers), and 



.8 


inclusion of higher order terms extends the validity into lower 
Reynolds numbers. 

Of course the basic premise of such a perturbation scheme 
has limitations. One such limitation applies Then, the neglecting 
of one of the terms changes the character of the governing equations 
drastically. Such a condition is called singular behaviour and 
examples of this include the boundary layers on blunt bodies which 
give rise to 1' .Alembert's paradox. 

To be able to apply perturbation methods one must he able 
to compare the orders of terms. It is helpful to introduce a formal 
concept of ordering. 

Ordering : 

Consider a class of function f(e) of a small parameter e , 
such that 

i) f is real and continuous 

ii) for each f there exists an interval for which the function 
is positive (of the type 0 <_ e <_ e Q ). 

iii) We take lim f(e) = 0 that is limit of the function 

e -*• 0 

as e approaches zero is zero. 

The order of these functions refers to the relative rate at 
which these approach zero. To be able to express these rates 
qualitatively, one compares the limiting behaviour of the functions 
f with a set of gauge functions. These gauge functions g(e) consist 
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of an ordered set of functions of e which are so f amiliar that 
their limiting behaviour can be regarded as known. 

If we take limit of f/g as e tends to zero there are three 
possibilities 



f = 0 [g] (2.5) 

If case (i) holds we say that f is of higher order (tends to zero 
faster) than g and write 

f < g or (2.6) 

f = o [g] (2.7) 

The set of gauge functions g Q ( £ ) is ordered such that 

g 1 ( e ) > gg(e) > gj(e ) ... . (2.8) 


One realizes that there are infinitely many sets of gauge 
functions possible. The choice of an appropriate set is quite a 
bit a matter of experience and judgement. Before we go further we 
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introduce two additional concepts those of asymptotic sequence and 
series. 

Asymptotic Sequence ard Series ' 


An asymptotic sequence is one in which the terms are arranged 

in an ascending order. Thus a sequence "f e ) ^ is asymptotic if 

$ , „ < 4> , that is 
n+i n T 


4> (e) 

-ffef-0 

e 0 


( 2 . 9 ) 


A series £ a^ $ (e) is termed an asymptotic series if 


{$ } is an asymptotic sequence. 


An asymptotic series £ a^ is an asymptotic expansion of 
a function f(x,e) to IT terms if ard only if 
IT 


f(x,e) = I a n (x) * n (e) + 0 [^ ] 


( 2 . 10 ) 

that is the remainder after IT terms tends to zero faster than <J> 


IT 


as e -*■ 0, This can also be put as 
TT-1 

f(x,e) = I ajx) 4> n ( e ) + 0 [4^] 


( 2 . 11 ) 


Here a *s are given by 

m ° J m— 1 


a = lxm 
m 

e ■+ 0 


f(x,e) - l a 6 
1 


( 2 . 12 ) 


Therefore 


a = !ia F 
1 e 0 i 


( 2 . 13 ) 
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a = lim 

^ e o 


a 6 
1 V 1 


( 2 . 14 ) 


The Perturbation Technique 

The perturbation technique consists of recognizing the parameter 
6 for whose limiting value the solution is known or can easily be 
obtained. The actual solution f(x,e) can then be thou gh t of as a 
perturbation on this basic (or first order) solution. The dependent 
variables are expanded as an asymptotic series whose first term is 
the basic solution and higher terms can be written as f^(x) 4> n ( e )» 
where ^4> n (e)} is an asymptotic sequezee. 

The next important step in the technique is to non-dimensionalize 
and normalize each variable by its characteristic value for the problem. 
This ensures that all the variables and derivatives in the equation 
are of older unity, and the relative orders of the various terms are 
determined by the gauge function which appear as the coefficients of 
the variables. By equating, terms of similar order we obtain a sequence 
of equations for functions f , that can be solved one by one. It is a 
property of this technique that only the basic equation can be non- 
linear, all the rest are necessarily linear. Once we solve for f n 
we can write the asymptotic expansion f & (x) 4>^(e). The radius of 
convergence of this series then determines "the domain of validity of 
this perturbation expansion. As happens sometimes, this domain of 
validity inc ludes the whole range of x we are interested in, and then 
the problem, is termed one cf uniform perturbation.. On the other hand, 
if the expansion obtained is rot uniformly valid, we term the problem 
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one in singular perturbation. The reason why the expansion breaks 
down in certain region of interest lies in the fact that the charac- 
teristic quantities used for noraali 2ation not proper fap ^ 

region and consequently the derivatives ard variables are not all of 
order one. This happens, for example, in case of boundary layer 
over a flat plate of length X, in which though the characteristic 
length for the 'potential flow' region is 1 , that for the 'bouMary- 
iayer' region is 5 , the boundary-layer thickness. This upsets the 
seperation of the equation into the sequence of equations for functions 
f n* T ° ° Tercome this difficulty we use different normalizing variables 
for this region and obtain an expansion valid for it. The procedure 
for the search of optimum nonnalizing quantities and the patching or 
matching of the two expansions at the interface of their regions of 

validity can be found in ary book on si^iar perturbation such as 
Cole ( 5 ). 

In our present study, because of uniform mechanisms of diffusion 

throughout the flow field, we have a uniform perturbation problem. More 
on this later. 

In the next chapter we pose the problem of a two-dimensional 
turbulent jet and develop the basic equation and boundary conditions 
for the problem* 
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Spread of a Turbulent Submerged Jet 
Velocity Variation along the of 

of a Submerged Jet 
Governing Equation and Boundary 
Conditions 



THE TWO-DBEHSIONAL TURBULEHT JET 


Introduction 

A free jet occurs v/he never a fluid is discharged from a 
nozzle or an orifice. For an inviscid flow a jet has one or more 
surfaces of discontinuities across which such parameters as flow 
velocities, temperature and mass concentration experience disconti- 
nuities. The presence of diffusion effects of viscosity, thermal 
conduction and mass transfer tends to dissipate the discontinuities 
so that a sort of boundary layer is formed. The instability associ- 
ated with the presence of points of inflexion in the velocity profile 
create eddies which move in a disorderly fashion and give rise to 
turbulence which further strengthens the momentum, heat and mass 
diffusion. It is experimentally determined that the jets are almost 
always turbulent even at quite moderate Reynolds numbers and even when 
a very smooth flow is maintained initially. We shall be concerned 
here in this chapter with constructing a model for a steady two- 
dimensional jet issuing out of a slit into stationary fluid (same 
fluid). We will also assume isothermal conditions so that thermal 
transfer and mass transfer effects do not enter the problem. 

Figure 1 shows simplified structure of a 2-D jet. Tie jet 
issuing out of the slit (gap * d r ) spreads as it moves in the x-direction 
entraining more and more fluid. Hie flew has some turbulence almost 
from the very beginning, though in the initial region there appears 
a potential core, in which velocity is substantially constant. This 





15 


potential core is eaten up by the viscous action and the turbulence as 
the flow proceeds downstream. Ibis region is known, in literature as 
the mixing region. Next comes the adjustment region beyond x = 4d 
in which the velocity profile adjusts and becomes self -similar in the 
main region. In this main region which starts at around x = 8d, the 
jet becomes similar in appearance to flow from a line source of infini- 
tely small thickness. The velocity profiles are self-similar to at- 
least the first order. This self-similarity arises from the fact that 
far away from the slit there appears to be no characteristic length 
dimension in the problem, except for the eddy size, which should not 
be relevant at the macroscopic scale. 

As the surrounding pressure impresses itself upon the jet 
stream from all sides, the pressure gradient is very small, and atleast 
to first order it can be taken as zero, as will be shown# The total 
momentum in the x-direction therefore remains constant to the first order. 

'Port hm arm studied the velocity profiles of an air jet emerging 
from a slot 0.03 m high with a velocity of 35 m/sec at distances x = 0, 
0.2 , 0.35 , 0.5 , 0.625 aid 0.75. The results are shown in Pig. 2. They 
show a continuous broadening of the velocity profile of the jet. The 
profiles become lower and wider. When plotted in dimensionless form, 
that is u/u vs y/y where u is centre-line velocity and y 
distance where velocity is half of u q , then for the main region they 
fall, more or less on the same curve. See Pig. 3# 





FIG. 3 _ DIMENSIONLESS VELOCITY PROFILE IN PLANE JET 




u = u + u 1 and 

(3.1) 

V = V + V* 

Irandtl (14) assumed that the rate of spread of jet (measured by the 
material rate of increase of jet half-width b with time) is directly 
proportional to the transverse component of fluctuating velocity, 
fhus 

~*v» (3.2) 

dt 

She similarity of the profiles across the jet also implies that the 
non-dimensional mixing lengths at the corresponding points should not 
vary with the axial location. The similarity of velocity profiles 
is obtained by scaling the transverse length by b, the jet half— width 
and the velocity by u Q , the centre line velocity (or the maximum 
velocity). Hie similarity therefore gives 

m = aod. 


(3.3) 
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where c is a constant of proportionality. 

The above result with the fact of similarity of velocity 
profile gives that the constant velocity lines are straight rays 
issuing out of the point at x = 0. If we take the velocity profiles 
in the fully developed region and extend the rays backwards, they 
converge to a point slightly upstream of the slit. See Fig. 4. 
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FIG. 4 .EQUAL VELOCITY LINES ( u/u m = CONST ) IN MAIN 
REGION OF SUBMERGED JET . 
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.Mpcity variation along the axis of a sutoerged ; 

Pressure gradient in a jet is zero atleast to the first order, 
hence the total momentum of fluid mass per second must be same in all 
cross sections of the j et. Therefore 

iA = Const. (3.12) 


/ P u 2 


Area 

Por a plane submerged jet we can write this as 

0 V f (t - ) 2 (3.13) 

O 0 X * 

Because of similarity of velocity profiles we have 


r = *(y/*) 

0 

b/x 

Hence p f (— ) 2 d© = 

J o V x 

(3-13 ) and (3.15) give 


Const. 


u x = Const. , 
o 7 


and therefore 


u = 
0 


Const. 


& 


Governing Equation and boundary conditions ; 


(3- 14) 
(3.15) 


( 3 . 16 ) 


Yorticity transport equation for mean-steady, incompressible, 
two-dimensional turbulent flow can be obtained from momentum equation 
as 

U(V 2 if>),x + v(v 2 i}>),y = |v\ + r_,yy + T_,xy - t w »xx - 


xy 


XX 




yy 

. ( 3 - 17 ) 
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where , 


U = v/U e , 

(3.18a) 

t = v/n e 

(3.18b) 

1 ! 

> 

( 3 . 19a) 

y s y*/d 

(3.19b) 


The Reynolds number 'R' is defined as 
P U d 

R * f (3.20) 


U e = Jet exhaust velocity 
d = Width of jet slit' 


x' and y' are physical coordinates, x’ being along the jet axis and 
y f normal to it. u and v are mean velocities along these axes. Prom 
the definition of stream function which integrates the continuity 
equation, 



( 3 . 21 ) 


V = 


djp_ 

dx 


( 3 . 22 ) 


x 1 s are the Reynolds stress components defined as 


t = - u' 

XX 


T = — U T V* 

xy 


T = - 

yy 


(3.23a) 

(3.23^) 


where 


(3.23c) 
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u T ' ai ’ e non-dimensional, fluctuating components of the velocities 
along the j ad y directions. The bar denotes time mean.' 

In the vortic ity equation, commadenotes differentiation# 

For a plane-jet the relevant boundary conditions are : 

At y = 0 


u = u 

0 

(3.24a) 

7 = 0 

(3.24b) 

sia 

1! 

O 

(3.24c) 

u ^ 0 

(3.25a) 

u» -> 0 

(3.25b) 

v 1 -* 0 

(3.25c) 

£U 0 

dy 

(3*25d) 


In the next chapter we will set up a perturbation scheme to 
solve equation (3.17) with boundary conditions (3.24) and (3.25) and 
obtain solutions upto second order. 
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The Perturbation Expansion 
Solution of First Order Equation 
Solution of Second Order Equation 



PERTUBBAIIOU EXPANSIONS OP SHE TWO-DMMSIOIIAI, JET EQUATION 


The Perturbation Expansion : 


The governing equation with boundary conditions was obtained 
in the last chapter as : 


U(V 2 ifO,x + v(v%),y = -J ip + t ,yy + t ,xy - t ,xx - t ,xy 

■tt 3y xx xy yy 


u = 

, and 

3y ’ 

(4.2a) 

Y = 

dip 
” 3x 

(4.2b) 

The relevant 

boundary conditions are: 


At y = 0 : 

U = U 

0 

(4.3a) 


7 = 0 

(4.3b) 


& = o 
dy 

(4.3c) 

As y -*■ co 

U-v 0 

(4.4a) 


u’-s- 0 

(4.4b) 


v' -¥ 0 

(4.4c ) 


+ 0 
dy 

(4.4d) 


In free turbulence the distribution of turbulent energy and 
dissipation are nearly uniform in the main body of the flow, and .unlike 
wall turbulence, are not dominated by gradient diffusion. A large 
scale convection takes place which transfers turbulence in bulk (25 )• 
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In the case of wall shear flow where the vortieity and the 
turoulence energy are generated at the wall and are then dissipated 
into the main flow by the viscous and turbulent action, there are 
two distinct regions of flow, one near the ,>cil in which viscous 
action predominates, and the other far sway from the wall where 
turbulence dissipation governs the flow. In such a situation we can 
not construct a uniformly valid perturbation scheme and the problem 
is necessarily one of singular perturbations. In this case of free 
shear flow, however, since the viscous and turbulent dissipation 
are nearly uniformly distributed, we may expect to construct a 
perturbation model which will apply throughout the flow region, dius 
the problem appears to be a regular perturbation problem. We will 
infact see that a regular perturbation scheme works, and in the 
following we develop such a schane. 


In this problan we scale the length x such that the new 
independent length variables X and Y become of the same order. We 
will use Reynolds number as the parameter, and expand the solutions 


about — = 0. 

ii 


X = A(r) . x 


(4.5) 


a(r) is a small parameter 


Y = y (4.6) 

fflae expansions of dependent variables can be written as 

iKx,y;R) = E^R) ^(x,Y) + E 2 (r) t|> 2 ( x > y ) + 0 [E 2 (r)] (4.7) 

U(x,y;R) = E t (R) * 1>y (x,Y) + e 2 (e) *2,Y (X ’ Y) + ° [E 2 (r) 3 (4 ' 8) 
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where E^(r) and E^r) are gauge functions belonging to an asymptotic 
sequence (E^Cr)} such that 

T E ^1 U) n 
lam -p TrY *0 

R -*- » V Rj 

Y/e next expand the Reynolds stresses in asymptotic series. For 
this purpose, we adopt the closure hypothesis proposed by Prandtl, 
which is known as Prandtl ! s new theory of free turbulence* For two- 
dime nsional flow, he approximates the Reynolds stress components by 


T xx 


t = 


§li. «i = 2E fi 

dy dx dx 


2 dU dU 


dU 


xy 


yy 


£ dy * dy K dy 

„2 dE dV dV 

dy dy dy 


(4.9) 

(4.10) 

(4.11) 


In these expressions by Prandtl 's new theory the eddy viscosity 
coefficient K - 5,2 ~ , can be considered independent of the y-coordinate, 
that is at any station x, it remains constant across the entire cross 
section of the flow. Prandtl obtained an expression for the eddy 
viscosity coefficient K by considering the similarity of non-dimensional 
velocity profiles at various cross sections of the flow which permits 
us to write 


d(TJ/U ) 


(4.12a) 


' d( ' y/b"J = Inv ’ ariarrfc f ' or a free 5 e ' t 

iO.sc the non-dimensional mixing leqgth (expressed in fractions 
of jet width) is same for corresponding points at different cross 


sections of flow 
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Therefore — = Invariant 
b 


(4.12b) 


Using (4.12a) and (4.12b) we can write for the eddy viscosity coefficient 
K for a free submerged jet : 


K = 


dU 

d y 


= b 


U 

o 



d(u/u Q ) 

d (y/b) 


or 


E = x 1 b U q where (4.12) 

X.^ is a constant. It has been experimentally determined that K 
remains approximately constant near the centre of the flow field but 
decreases somewhat towards the edges. Thus the results obtained 
with this hypothesis match the experimental results near the centre 
but show some discrepancy towards the edges. 

In chapter III we have shown that the centre line velocity 
decreases as l/ and the j et half -width increases as x. Thus 

K = X, U b 
1 o 

= G C /X (4. 13 ) 

1 

where G^ is an appropriate gauge function undetermined so far and 0 
is a constant. 

Using the expansions (4.7) and (4.8) and definitions (4.5) and 
( 4 . 6 ) of independent variables* we obtain the various terms of the 
governing equation (4.1 ) as 
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U(v 2 4i),x + V(V^),y * e£vU. * ' , \ 

1 '*1,Y V 1,m *1,X*i,m' 

S 1 V U 1,Y *2,XYY + ^2,Y " %,X *2,YYY 


+ E. 


“ ^2,X ^1,YYY^ + ° ®- E 1 E 2 & ^ 


\ (7 4 » 


E 


__L 

R " r 1,YYYY 


E 


^I.YYYV 3 + ° EjT 3 


xy 


XX 


yy 


” & i 0 ^ (e i *i,ir + E 2 ^ 2 ,yy } + 0 [& i e 2 3 

= Vi A (2C /l * 1jXy ) + o [E 1 A ] 

= B 1 & 1 A(2C Jx ^ 1yXy ) + o [E 1 & 1 A] 


(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 


On differentiation of (4.16), (4,17) and ( 4 . 18) we get 


2y»yy 


xx,xy 






= E 1 & 1 ( 0 ^^ 1>YYYY ) + ^\(C -/X ^ 2>YYYY ) + 0 [ff^J (4.19) 
2 

“ E i e i 4 { alaf (2 ° 5 *i,rr )! + 0 [E i 8 i ^ (4 - 20 > 

2 

= E i 8 i 4 { alf +0 ! E i 8 i fi2 i 


= E 


2 r d 2 


>lV f ^ (C ^^1,YY )}+0[E 1 & 1 A 3 


(4.21) 

(4.22) 


In free turbulent' flows the turbulent stresses are much larger 
than viscous stresses in the entire region of flow, that is, the viscous 
stress terns are of- a- higher order than the Reynolds stress terms (20). 
fflie orflly important contribution from the Reynolds stress terms for the 
lowest order equation comes from the lowest order term in the 

rest being of a higher order, that is, being much smaller in magnitude. 
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Hie lowest order equation must therefore contain the inertia terms 
and the Reynolds stress terms t Thus the orders of the lowest 

order terms in the expansions of the two must be same, and we obtain 

* E 1 G 1 

We note that 1 to satisfy the condition on the jet axis, 

that is, U U Q as Y 0. Without loss of generality, we can set 


E.j = 1 , and thus obtain (4.23) 

G^ ^ A . We can again, without loss of generality, set 

& 1 = A (4.24) 

The lowest order equation thus reads : 


t 'P 
1 ,Y 1,XYY 



4 > 

1,YYY 


= C ^ 


1,YYYY 


(4.25) 


We see that even to second order the only contribution from 
Reynolds stress terms comes from t , since the order of the 

xy>yy 

second term in expansion of t (4.19) is lower than the lowest 

order terms in the other Reynolds stresses. We take the viscous stress 
terms to be of second order. The second order equation must therefore 
contain the second order inertia terms (4.14), the second term in 
the expansion of f (4.19) and the lowest order term in the 

xy,yy 

expansion of viscous stresses (4.15). Thus the orders of these terms 


must be the same and we obtain 

E 1 

E 1 V * R~ * E 2 G 1 


(4.26a) 
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Por the sequence of gauge functions E^, Eg we will choose 
a power series in the small parameter A which gives, Eg ' u - A. Ihis 
choice can be justified by the argument that if the Reynolds stresses, 
other than r xy ^ are to contribute even to the third order equation 
they must be of the same order as the third order terms in the 
expansion of inertia forces. For this we must have Eg a, a. Without 
loss of generality we can put, 

Eg = A (4.26b) 

We obtain a relation between R and A from (4.24) and (4.26)- 
a and b 


A a, — , and without loss of generality 

/R 


A = 


/R 


(4.27) 


Hence the first and second order equations from (4. 13 )» (4.14), 
(4*15) and (4*19) turn out to be 


1,Y 

^1,XYY 

* +i,i 

^1 ,YYT 

c Sx 

*1 

,YYYY 

r,T 

^2,XYY 

+ ii 2jY 


h,x 

^2 

,YYY ‘ 


* + 0 ’ / X 4> 

*1,YYXY 2,IYYY 


(4*28) 


(4*29) 


Hie corresponding boundary conditions are : 


= U 

(i) 

1-,Y o 


= 0 

(ix) 

1,x 


1 =0 

(iii) 

1,YY 



At Y = 0 


(4.28a) 
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As 1 -*■ 00 

V*° 

(i) 



^IjYY'*' 0 

(ii) 

(4.28b) 

and for conditions on we have 



At Y = 0 

it =0 

V 2,X 

(i) 



*2,n “ 0 

(ii) 

(4.29a) 

As Y 00 

*2,Y - ° 

(i) 



tl i -*■ 0 

^2,YY 

(ii) 

(4.29b) 

Solution of 

First Order Equation : 




The equation we obtain on integrating the first order equation 
(4,28) once with respect to Y is of the nature of boundary layer equation 
with zero pressure gradient and. reads as 




ip 

1,Y %XY 








1 ,YYY 


(4.30) 


It should be noted here that since the vorticity equation is 
obtained by differentiating the momentum equation where in the pressure 
gradient terms disappear, on. partial integration of the vorticity 
transport equation, we should have a function of integration which 
should depend upon the pressure gradient. Since the pressure gradient 
along the j et axis is zero to the first order the constant of integration 
is also zero. 

This is the same equation as obtained and solved by Soertler for 
tile two-dimensional turbulent jet based on Prandtl T s second hypothesis (8) 
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with the difference that this is a non-dimensionalized and scaled 
equation. Adopting the same approach we introduce h and f ^ (n), such 
that 

^ = A ’'Xf^Cn) , where (4.31a) 

11 = T (4i3lb) 

A and B are constants. 

Using (3.3 l)a and b we get for the various terms of the first 
order equation (4.30 ) 


* f'(n) 

1 »* € 1 

_2 

(4*32a) 

. f"(") 

X v 5 T 1 

2 

(4.32b) 

+1.YXY ^*1 (n > 

(4.32c) 


(4.32d) 


(4.32e) 


Substituting (4.32) in (4.30) we obtain 

_ ti- f ; 2 (n) i£- „ f'(n) f"(n) - ^ f”(i) *,(") 

2 X 1 r 1 1 2X 

x 2 1 1 X 2 

On simplifying and rearranging we get 

f» T 2 (n) + f;(i) f,(n) + “ = 0 


(4.33) 
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As we have two free constants A and B we can for convenience 


put 


A = 4BC 

with wh ich the first order equation reduces to 
2f^ 2 (n) + 2f»U) f^n) + f»»(n) = o 


(4.34) 


(4.35) 


This equation must be solved subject to the boundary conditions 
(4*28) a and b which reduce to the following 


n = 0 


f ;- 1 


f i = ° 


f ' 0 

1 


(4.36a) 

(4.36b) 

(4.37) 


Integrating (4.35) once we obtain 


2f t (n) f'(n) + f»(n) = Const. = 0 


Integrating this once more we obtain 


2 

f^(n) + f *(n) = C^n + C^ . Using boundary conditions 


(4.36) and (4.37) we obtain 


c i =° 


°2 = 1 


Therefore, f 2 (n) + fj(n) = 1 (4.38) 

This equation is same as the one obtained for two-dimensional 
laminar jets and solution to this is 

C t.i- 


j.R 

> 4 RY 


& 47747 

Acc. *v>. *.*•**“* 
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£ ^ (n) = tanh t) 


and i|i = A /E tanh p 

Hence we get for the velocities 


U 


i - V = i (l " tanh ^ n) 


T = - 


*1,X 


AA 

2^T 


(4.39a) 

(4.39b) 

(4.40) 


[tarh p — 2 n ( 1 - tanh 2 p ) 3 (4.41 ) 


Ihus we see that the non-dimensional veLoeity profile is self- 
similar to the first order atleast. 


For evaluation of the constants we use the condition that 
pressure gradient is zero to the first order and therefore the first 
order momentum flux in x direction must remain constant. 


Momentum flux 


= J = 



dy 


+co 


Therefore 


2 ^ 
d tJ / (ip )* dy to the first order 

6 I j I 


= P d U 2 / (l - taifc 2 p) 2 dy 

G A. 

mm 00 

o p 00 2 2 

= P d U A1 I (l - tanh n ) dn 

^ mm CO 

let £ = tanh n 

2 

Therefore d£ = (l - tanh p) dt 

+1 

Hence J = P d U 2 A^b / (i - £ 2 ) d£ 

= P d U 2 A 2 b(4/3) 

6 
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Therefore — A^B = ~ = 1, and 


P d tJ 


A = 


/W 


(4.42) 


Prom (4.34) and (4.42) 

k _ € 


C 


4B 8B^B 


(4.43) 


Ebw the only unknown constant is 3 and this is evaluated 
expe rime nt ally. 

To investigate for the contribution of pressure gradient term, to 
the second order equation let us consider the y-mamentum equation 

.2 Tr 3 t 9t 


T j^_ + ? 32: = _m + i + £ 1 ) + + __zz, 

u <V«. V ' r\TT T> ^ O ‘ nJ ' rv w ^ rw,. 


3Y 

ax ■ ay 


9y ' E Vx 2 ' 9y 2 ' ' 3X 9y 


(4.44) 


If we use the expansions of U., Y etc. in this equation the lowest 
order equation becomes 


- ^ v v- - 4>. , ™ - 2C ^ wv) = “ 


1,Y T 1,XX r 1,X 1 ,XY 


1,XZY' 3Y 


(4.45) 


Thus ^A 2 


Since Y is of order one we can write 

2 


P * A 


We next find 


3p _ 3£ 3Y 
3n 3 Y # 3n 

_X ££ 

" B 9Y 
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p is a function of n and x therefore 


|E= A^B 

. ax ax 


- & (IB. la + IB) 

v 9n sx ax' 

= a (IB _ H IB) 

L v ax b ay' 


The terns in the bracket are each of order h so that 

IB n, &3 
ax 

Thus we see that is of a order higher than second order 

ax 

also, or is smaller in magnitude and can be neglected from the second 
order equation also. 

Solution of Second Order Equation ; 

Integrating the second order equation (4.29) once with respect 
to Y we obtain 

^1,Y ^2,XY + ^2,Y ^1,XY ** *1,X ^2,YY “ ^2,X ^1,YY 


t „„. r + c 

V-< 


,IYY 


2,YYT 


(4.50) 


Since the pressure gradient term has been shown to be of an 
order higher than the second order, we can take the function of 
integration to be zero. 

The first order terms are given by 

ih = a /X t'arih p with p = BY/X (from eqn. 4.39) 
1 
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V = f 0-W.o 

2 

_ 2ab 


* 


1 ,IY ~ tanh n ( 1 - tarih 2 n ) 

3 


2AB' 


1,YYY j2 n (l - 'tarn 2 n) - (l - tarih 2 n)] 


= _A_ 
1?x 2/r 


[-tarib. n - 2 n (i - tanh 2 n )] 


(4.51 ) 


AB 


^» xy = 2X^ [4n tanil n (l " taJah n) - (l - tarih 2 n) 3 


AB , 2 

PP== (1 - tarih n) (4n tarih n - i) 


Substituting from. (4.51 ) in (4.50 ) we obtain 


■— (l - tarih 2 n) * + (4ri tarih n - l) (l - tanh 2 n) ip 

* 2XvY 


2,Y 


A 


(tarih n ^ 2n (i - tanh 2 n)] * 0 w + ■==- tarih n (l-tanh n) ip v 

2/? ~ 


C 




2 ,YYY x 2^ 


2_ ^ ( 1 - tarih 2 ri ) (3 tarih 2 n - 1 ) 


(4.52) 


Multiplying through out by 


x 2 /x 

r we get 

AB(l - tanh n) 

tarih ri 


- ri 3 


*2,XY + X ( 2ritaj± * %j ~ B" ^2 (l-tarih 2 n ) MJr2 ’ YY 


ci 3 ^ 

2,X AB(l-tarih 2 n) 


+ 2B tarih n _ 2? ^F~ — = 2B 2 (3 tanh 2 ri - 1 ) 


(4.53) 


She solution of this partial differential equation would bee ctae 
simpler if we could reduce it to an ordinary differential equation 
by similarity arguments. To explore for similarity in the second order 
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equation we let 




2 


n 


2 


A X a f (n ), and 
BY 


(4.54a) 

(4.54b) 


Using (4.54) we get for the various terms of the second order 
equation (4.53) 


\p 

2,YY 

^2,YYY 

^2,X 

^2 } XY 


‘ AB x 3 ' 43 f^(n 2 ) 

= ab 2 x a-2b f"(n 2 ) 

= AB 3 X^ 3b f| ,, (r» 2 ) 

= a A X 3 " 1 f^) " A Xa " 1 \ f 2^2' > 
t= a A B X 3 " 13 " 1 f J(ti 2 ) - AB X^’ 1 n 2 

- ab x a “ b_1 f^(n 2 ) 


f 2 (n 2> 


(4.55a) 

(4.55b) 

(4.55c) 

(4.554) 

(4.55e) 


Etot (4,55) ani (4.53) l» see that for the coefficients of 

* • * ; * (*„ T am 4, ™ to he independent of X we 

2,XY’ 2,Y’ 2,XY 2,X 2,YYY 


must have 


a «x b ■ - 1 = -2 
a - b * -1 

a - 2b = -2 

a - 1 = -1 

a - 3b = -3 


(4.56a) 

(4.56b) 

(4.56c) 

(4.56a) 

(4.56e) 
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These conditions are satisfied if we pu t 


a = 0 


(4.57a) 

b = 1 


(4.57b) 

so that 



n 2 = n 


(4.58) 

and we obtain 



4*2 = A f 2 (n) where 

„ by 

n = r 

(4.59) 

With these expressions for 

and n we get 


4> 0 v - f ’( 0 ) 

2 j x X 2 


(4.60a) 

^2,YY = ~ *2^ 


(4.60b) 

*2,YYY =|r f 2 ,,(ri) 


(4.60c) 

^2,X = ~ x” *2^ 


(4.60d) 

j, ss MU f’(n) - — 

*2,T£ x 2 V ' x 2 


(4.60e) 

Substituting from (4.60) in (4.53) we get the 

second order 


equation as 

-BA [nf”(n) + f^(n) 3 + A B f^(n) (2n tanh n-|) 

r tanh r) __ f «( n ) - 2ABn tanh n . f^( n) 

2(1 - tanh 2 n) 

r B 2 0 f „,( n ) = 2B 2 (3 tarh 2 n - 0 

~ l r ^ , 2 v J 2 v ' 

(1-tanh n; 
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% 

AB taah n _ f „( n ) - - AB f 'K 

^ ■ p ^2 ^ 7 2 t 

2(1 — tanh. n ) 


^ *i •* t o-hh. ^ ^ 


f ni(n) 


o 2 * 

= 2B (3 tanh n “ ' 


To simplify this equation, multiply 


y throughout by 


. ,2 _\ 4(1 - tarfa JOj 

1 — tarih jn ) or by — ■ ■ 43 

A * 


as G = 4 B 


This gives 


, w fi (, - ta* 2 n) *aW 
f»»(q) + 2 tanh n fg^' + 


+ JSL (i - tanh n) (3 t arl1 


2 = 0 

n - 1 / 


(4.61) 


We are 


interested in velocity Ug 


*icfc is proportional to f 2 ‘ 


Thus we can reduce this 
second order by putting 


third order differ® 


atial equation to one of 


fi(ri) = z(n) 
2 


(4.62) 


■This gives, 


2 ii(b) + 2 tanh b Z'(h) + 


2 n ) z(n) 


8 B (1 _ tanh 2 h) (3 talk 2 n - l} 
A 


(4.63) 


At n = 0we have 


z „(0) + 6Z(0) -I 




(4.64) 


As b -*■ 00 we have 


(4.65) 


z ,,(«) + 2Z'{ m ) ~ 0 
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Let 


2 tanh n = g^q) 


6(l - tanh n) = g 2 (n) 


8B 


/ 2 p 

(1 - tanh n) (3 tanh n - l) =^(n) 


(4.66a) 

(4.66b) 

(4.66c) 


Ehen the second order equation reduces to, 

Z"(n) + g^n) z'(n) + g 2 (n) z(n) + g^(n) = 0 

The boundary conditions to be satisfied are 
At n = 0 : Z* = 0 

As n -»■ 00 : Z -v 0 


(4.67) 

K68) 

(4.69) 


Obtaining an analytical solution to this equation is a difficult 
task because of the variable coefficients g^(n), g 2 (n) etc. involved. 
We will resort to numerical techniques using the principle of super- 
position, outlined below. 

The Numerical Technique 

An analytical solution of a second order equation of the type of 
(4.67) will be of the form 

Z - °1 %, + °2 + h (4 - 70) 

where an d Z^ .-.are solutions of the homogeneous part of the 

differential equation and Z p is the particular solution. 0 1 and G 2 
are constants. Because of the linearity of the equation we can solve 
for , Z H2 ani Zp independently and then add the results, after 
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determination of the constants and to satisfy the boundary- 
conditions. 

The homogeneous part of the differential equation (4.67) reads 


Zg(n) + g.,00 Zg.(b) + g 2 (n) ZgCn) 


(4.71 ) 


For the initial values we will organize this as 


V 0) 


(4.72a; 


(4.72b) 


Z H2 (0) " 


(4.73a) 


(4.73b) 


For the particular solution the entire equation has to be solved, 


that is, 


z”(n ) + g^n) Z^(n) + g 2 (n) Zp(n) + ^(n) 


(4.74) 


The initial values for this will be taken as 


Zp(0) = 


(4.75a) 


(4.75b) 
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Hie constants and in equation (4.70) have to be evaluated 
from the overall boundary conditions as follows 


At n = 0, if we have Z = a we can write 

Z(0) = a = 0 1 Zg 1 (o) + C 2 Z H2 ( 0 ) + Zp(0) 



(Therefore = a 

As h + ® , if we have Z = b we can write as 

z(co) = b = c 1 + c 2 z H 2 (») + Zp(») 

t> - C 1 2L (<») - Z_(«) 

(Therefore C 2 = Z~T»J 


(4.76) 


(4.77) 


For our problem the value of Z as q -*• «° is zero, that is, 
b = 0. (The initial value z(o) is not known, however from equation 
(4.64) we see that to start the integration procedure we can take 


z(o) 


8B 

6A 


It was also found that with the initial values as given by ( 4 .73), 
the solution of the homogeneous equation was a trivial one, that is, 
the value remains zero throughout. (This caused the solution Of C to 
blow up. To avoid this Z ff2 (o) was given a small value, of the order 
of 1 x io“ 5 . 2 ^ , Zjj 2 and Z^ have been solved by the Eunge Kutta 

method. Hie results obtained have been discussed in the next chapter. 



CHAPTER. 5 


RESULTS Alffi DISCUSSION 




RESULTS MD hlSOJSSIOIT 


In the last chapter's we developed the equations for first and 

second-order velocity distribution using perturbation theory arguments 

One of the main points that emerged in the analysis is that viscosity 

effects which are negligible compared to the Reynolds stress effects 

in the first order become significant at the second order. Also it 

was shown that even to second order the pressure variations are 

negligible compared to the effects of the other terms in the momentum 

2 4 

balance equation. Since for the range 10 < Re < 5 x 10 the second 

order correction is small, the effect of pressure variation on -the 

2 

jet in this region is really very small. Only for Re < 10 ,, when the 

second order correction is substantial and indicates that the third 
term may need to be considered, can the pressure variation have any 
possible significance. 

The first order velocity profile, as was expected, is independent 
of the Reynolds number and of x-station. 'This infact represents 
the solution in the limit x -»- » and Re <» . Eor moderately high 
Reynolds number and for moderately large x/d (as opposed to large 
x/d and Re) this profile needs correction. We calculated, using 
the procedure outlined in the previous chapter the second order 
correction to the basic velocity profile for four x-stations ( 10 , 15, 
20 and 25 jet widths from the slit) and five jet-Reynolds numbers 
(lO 2 , 10 3 , 5 x 10 3 , 10 4 ani 5 x 10 4 ). The results obtained are 
presented as the following figures. 
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Figure 5 shows the second-order corrections across the jet for 

Re = 5 * IX) 4 and x/d = 25. It should be noted that the effect of 

viscosity is to reduce the velocity near the centre line and to 

increase the velocity at the edges. This is. the typical effect expected 

of viscosity-tc diffuse the velocity gradient so as to make the 

velocity differentials •smaller. Figures 6 and 7 show the second 

3 

order corrections across the jet for Re = 10 and x/d = 25; and 
3 

Re = 10 and x/d = 15 respectively. 

Figure 8 shows the variation of corrected centre-line velocity 
with Reynolds number at a given x-station (x/d = 20 ). The effect 
obtained is anall, increasing with decreasing Reynolds number. As 
Reynolds number increases the corrected centre-line velocity tends 
to its first order value. This is also in line with the explanations 
since the decreasicg Reynolds number mean the increasing viscous effects. 

Figure 9 shows the variation in the second order correction to 
the centre-line velocity with x-station for a typical Reynolds number 

’T 

(iCT )* ELgure 10 also shows the percentage correction in the normalized 

velocity (centre-line velocity taken as unity) at a given non-dimensional 

4 

y-location (n = 1 ) and a given Reynolds number =10'. Both the figures 
exhibit a decrease in correction as x increases. It can be inferred 
that the perturbation results obtained are valid for x/d between 8 
and 25. Beyond 25, the unperturbed results are good enough, while 
below 8, the perturbation scheme breaks down because of the absence 
of full development# 
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FIG. 5. VARIATION IN CORRECTION DUE TO VISCOSI1 

REYNOLDS NO: BxlO 4 STATION (x/d) . 25 






FIG. 8- VARIATION IN CENTRE 








ABSOLUTE CORRECTION 
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Figure 11 plots the corrections in the velocity profile, noimalized 
with respect to the corrected centre-line velocity. This is plotted 
for a given x-lo cation (s/d = 20 ) and a given Reynolds numb er (5 * 10 5 ). 
The significance of this curve is the fact that this brings out the 
nature of the second order corrections as would be observed in 
experimental results where the curves are plotted with the velocities 
normalized with respect to the observed centre-line velocity. 

Comparing figure 5 with the data reported in literature (20 ) we 
note that it does infact represent the experimental results for 
values of non-dimensional n = y/d less that 1.2 (n ^ Y/X less than 
approximately 1*15). For higher values of n the result breaks down. 

This is due to the difficulty associated with the closure hypothesis. 

The closure hypothesis due to Prandtl which we use here asserts that 
K = £ ~ is constant across the width of the jet. This is not really 
true for n greater than approximately 1.15 beyond which the value 
reduces. Another restriction is imposed by the fact that the one 
■undetermined constant involved in the analysis is dependent on the 
jet size and is not a universal constant (lO). 

2 4 

These results are valid in the Reynolds number range 10 < Re < 10 

and for x— stations 8 < x/ d <25. For x/d > 25 and Re > 10 the first 
order velocity profile (based on infinitely large Reynolds number 
and f ull development assumptions) can be taken as valid. For 

s/d < 8 and Re < 10 2 , the jet model breaks down and the results obtained 

•,.l; ’ - ... •.•.•ir.ontfcl 

litre linger valid* io r*-~ 



The dots circles in the graphs do not indicate experimental 
values* 
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